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We report a first-principles study of (BaTi03) m /(BaO) n superlattices for a wide range of peri- 
odicities m/n. We show that such a system develops a polar zone-center instability for sufficiently 
large m/n ratio, which can be understood, at least qualitatively, from a simple electrostatic model 
and should lead to a ferroelectric ground-state. However, the analysis of the phonon dispersion 
curves also points out the appearance of stronger antiferroelectric instabilities at the zone bound- 
aries around m — 4, before the critical ratio for ferroelectricity is reached and which still domi- 
nate beyond it. The dominant character of the anti-ferroelectric instability is explained from the 
depolarizing field which hardens the ferroelectric mode. This analysis allows us to predict that, 
(BaTiC>3) m /(BaO)n superlattices should present an antiferroelectric ground state for m larger than 
4, which should smoothly evolve to a multidomain structure for increasing m values and only become 
ferroelectric for large m. 
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I. INTRODUCTION 

During the recent years, numerous works have been 
devoted to the study of interfacial effects in ferroelec- 
tric nanostructures. 1,2 Both theoretical and experimen- 
tal works have focused on different types of ferroelec- 
tric multilayers and superlattices including (i) ferroelec- 
tric ultrathin films between metalic electrodes (FE/Me 
interface), 3-7 (ii) superlattices combining two ferroelec- 
tric materials (FE/FE interface), 8-10 and (iii) superlat- 
tices in which a ferroelectric alternates with an incipient 
ferroelectric (FE/iFE interface). 11-15 These studies high- 
lighted the fact that three main factors govern the physics 
of multilayers: epitaxial strain, electrical boundary con- 
ditions and interfacial effects. 

The case of ferroelectric/insulator superlattices (FE/I 
interface) in which a ferroelectric alternates with a regu- 
lar insulator was only marginally addressed in the litera- 
ture. This might be due to the fact that no enhancement 
of the ferroelectric properties is expected in such systems 
due to the detrimental effect of the regular insulator on 
the ferroelectric properties. For some time, only the in- 
terface between AO and ABO3 oxides, where A is an 
alkaline-earth atom such as Sr or Ba, attracted some in- 
terest since AO oxides are playing a key role as buffer 
layer in the epitaxial growth of perovskites directly on 
a silicon substrate. 16 A first-principles study was previ- 
ously reported by some of us for (BaTi03) m / (BaO) n and 
(SrTi03) m /(SrO)„ superlattices grown on Si, but only 
the paraelectric periodicity n=6 and m=5 was consid- 
ered and the ferroelectric properties were not explicitly 
discussed. 

Nevertheless the focus on these particular FE/I inter- 
faces has lately increased. The Ruddlesden-Poper series 
can be considered as an example of this type. 18-20 This 



renewed interest is fueled by the possibility of modify- 
ing in a controlled way the delicate balance between all 
the instabilities present at the bulk ABO3 perovskites, 
that might be ferroelectric (FE) , antiferroelectric (AFE) , 
and/or antiferrodistortive (AFD) phase transitions by 
changing the composition and/or periodicity of the se- 
ries. The perspectives of stabilizing new exotic phases 
not present at the bulk parent compounds is a strong 
driving force for these studies. 1 

In this paper, we study the case of (BaTi03) m / (BaO)„ 
superlattices epitaxially grown on a SrTi03 substrate as 
a prototypical example of FE/I superlattices. We first 
consider the predictions of a simple electrostatic model 
as reported in Ref. 11. Then, we present first-principles 
calculations of the phonon band structure of the paraelec- 
tric (BaTi03) m /(BaO)„ superlattices for various layer 
thicknesses m and n, and discuss the effects of the pe- 
riodicity on the vibrational properties. We show that 
the predicted ground state differs from that of the sim- 
ple electrostatic model, and explain the reason for the 
discrepancy. 



II. TECHNICAL DETAILS 

First-principles simulations were performed in the 
framework of the Density Functional Theory (DFT) as 
implemented in the Abinit package. 21 All the results 
were calculated through the local density approximation 
(LDA) for the exchange-correlation energy. An energy 
cut-off of 45 Ha was used for the planewave expansion. 
The Teter parametrization 22 for pseudopotentials was 
employed where the following orbitals were treated as 
valence states: 5s, 5p and 6s for the Ba atom, 3s, 3p, 3d 
and 4s for the Ti atom and 2s and 2p for the O atom. 
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From the smallest to the biggest superlattices, tetrago- 
nal Monkhorst-Pack meshes from respectively 6 x 6 x 3 to 
6x6x1 were considered for the Brillouin zone sampling in 
order to produce accurate results. The vibrational prop- 
erties and Born effective charges were calculated using 
the Density functional perturbation theory (DFPT). 23 

BaTiC>3 can be epitaxially grown on BaO. 24 
The epitaxy is such that ABO 3 (001)||AO(001) and 
AB0 3 <110>||AO<100>, which means that the BaTi0 3 
unit cell is rotated 45° around the (001) BaO direction. 

The superlattices were simulated by means of a super- 
cell approximation with periodic boundary conditions, 
so short-circuit electrostatic boundary conditions across 
the whole supercell are naturally impossed. To define 
the cell that is periodically repeated in space, a generic 
formula (BaO-Ti02) m /(BaO) n was used and labeled by 
m/n, where n is the number of BaO oxide atomic layer 
and m the number of BaTi03 formula units. 1 '' 21 

The epitaxial strain was treated by fixing the cubic 
in-plane lattice constant of the superlattices. Only the 
case of an epitaxial strain corresponding to an hypothet- 
ical SrTiOa substrate was considered (theoretical LDA 
relaxed cubic cell parameter of SrTiOa: a = 3.84 A). For 
each thickness, we performed structural optimization of 
the atomic positions and the out-of-plane cell parameter 
of the superlattice in its highly symmetric phase corre- 
sponding to the tetragonal space group P4/mmm (No 
123). 

To establish the notation, we will call the plane parallel 
to the interface the (x, y) plane, whereas the perpendic- 
ular direction will be referred to as the z axis. 



III. ELECTROSTATIC MODEL 

As explained in Ref. 2,11, the energy of a ferroelec- 
tric superlattice between materials 1 and 2 under short- 
circuit boundary conditions can be estimated within a 
simple model which assumes homogeneous polarization 
in each layer and neglect interface corrections, as the sum 
of the internal energy for the two compounds forming the 
superlattice and an electrostatic correction, 

E(P U P 2 ) = mU x {P x )+nU 2 {P 2 )+E elec {m,n,P x ,P 2 ). 

In first approximation, the internal energies U\ and U 2 
can be determined at the bulk level in the absence of 
a macroscopic electric field, while the term E e i ec takes 
into account the electrostatic energy cost related to the 
polarization mismatch between the two types of layers. 
From now on, we will assume that material 1 is BaTi03 
(BTO) and material 2 is BaO (BO). 

A. Bulk internal energies 

In their relaxed cubic structures, we obtained the lat- 
tice constants obo=5-46 A and aBTO=3-94 A for BaO 



and BaTi0 3 respectively, which are in good agreement 
with previous LDA results. 17 In the NaCl structure, a 
cell parameter of 5.46 A in the BaO bulk corresponds 
to a Ba-Ba distance of 3.86 A. Therefore, when grown 
on a SrTi03 substrate, assuming a theoretical cubic in- 
plane lattice constant of 3.84 A, BaTi03 is under com- 
pressive strain and becomes tetragonal, as well as BaO 
(the epitaxial strain being applied on the Ba-Ba distance) 
in a smaller extent. 26 After relaxation with this epitaxial 
constraint, the out-of-plane cell parameters of the para- 
electric references are c BO =5.48 A and Cg TO =4. 01 A for 
BaO and BaTi03 respectively (where the superscript "0" 
refers to the fact that it is the c-parameter at zero polar- 
ization). The associated Born effective charge and elec- 
tronic dielectric tensors are reported in Table I. These 
values are comparable to those in the cubic structure, 
with a small anisotropy due to the tetragonal symmetry: 
most components of the Born effective charges and elec- 
tronic dielectric tensors are smaller in the z than in the 
x/y directions. The frequency of the lowest transverse 
optic (TO) mode are also reported in Table I together 
with their mode effective charges. In this constrained 
paraelectric configuration, BaO is stable but the mode 
polarized along the z direction has a slightly lower fre- 
quency than those along the x/y directions. 27 BaTi03 
shows a ferroelectric instability along the z direction with 
an imaginary frequency of 180i cm -1 and an associated 
mode effective charge of 8.47 e while in the x/y direc- 
tions, the frequency of the soft mode is shifted to real 
frequency (195 cm -1 ) and its mode effective charge is 
slightly reduced (8.00 e). 



BaO BaTi0 3 
x/y z x/y z 
~a 3.84 5.48 3.84 4T0T 

Z*(Ba) 2.83 2.78 2.81 2.77 
Z*(Oi) -2.83 -2.78 -2.27 -5.29 
Z*(Ti) - - 7.48 6.80 
Z*{0 2 / 3 ) - - -2.13/-5.90 -2.14 
e°° 4.39 4.38 6.65 6.45 

ujto 118 74 195 180i 

Z* 3.14 3.08 8.00 8.47 



TABLE I: Cell parameters (a, in A), diagonal components 
component of the Born effective charge tensors (Z*, in e), 
diagonal components component of the optical dielectric con- 
stant tensors (e°°), frequency of the lowest TO modes (ujto, 
in cm' 1 ) and associated mode effective charges (Z*, in e) for 
the bulk paraelectric BaO and BaTi03, under the constraint 
of an in-plane lattice constant imposed by an hypothetical 
SrTi0 3 substrate (3.84 A). Oi refers to the O in the BaO 
plane while O2/3 stands for the two equivalent oxygens in the 
Ti02 planes. The diagonal components of Z* and e°° along x 
and y are equal by symmetry. 

The compressive epitaxial strain imposed by the sub- 
strate favors a ferroelectric c-phase 28 and, therefore, 
BaTi03 effectively becomes an uniaxial ferroelectric. In 
what follows, we shall assume that all the polarization 
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and electric fields are directed along z. A full relaxation 
of both the atomic coordinates and lattice vectors gives 
rise to a ferroelectric ground state with a spontaneous po- 
larization of 38 /iC/cm 2 and a relaxed out-of-plane cell 
parameter of 4.07 A. The difference of energy between 
the ferroelectric ground state and the paraelectric refer- 
ence is A£'=16 meV which corresponds to the depth of 
the double well associated to the evolution of the poten- 
tial energy with the polarization (see red curve in Fig. 1). 
In principle, the determination of the shape of the double 
well potential energy, U(P), would require calculation at 
constrained P as proposed in Ref. 29. Here we used a 
more approximate but still reasonable approach. First, 
we froze in different fractions of the pattern of atomic dis- 
placements £btOj defined as the difference between the 
atomic positions in the relaxed ground states and the 
paraelectric reference structure. Second, for each frozen 
fraction of £btOj the c cell parameter was relaxed in or- 
der to include the strain effects in the internal energy. 
Then, the corresponding polarization was approximated 

by -Pbto = ^' ^n" ' Zi ■ wn ere ^bto.i is the displacement 
of atom i with respect to its high symmetric position, 
Z* its Born effective charge [with the geometry used in 
the present work, the (z, z) component of the Born ef- 
fective charge tensor] , and Vt the volume of the unit cell. 
The "0" in the superscript makes reference to the fact 
that this polarization is estimated assuming zero macro- 
scopic electric field. Finally, the internal energy is plotted 
as a function of the zero-field polarization (red curve in 
Fig. 1), and fitted to a polynomial expansion 



The only way to polarize BaO is within the subspace 
spanned by its TO modes. In the superlattice, the BaO 
layer is expected to be polarized by BaTiOa along the z 
direction and so the evolution of its associated internal 
energy [/bo was determined by freezing the pattern of 
displacements of the TO mode along the z direction given 
in Table I. As for BaTi03 the out-of-plane cell parameter 
was relaxed for each amplitude of the frozen pattern of 
displacements. Since BaO is not ferroelectric, to polarize 
it has an energy cost and the shape of the internal energy 
corresponds to a single well (blue curve in Fig. 1) that 
can be approximated as 

Ubo(Pbo) = ^bo-Pbo + ^bo-Pbo' ( 3 ) 

where Abo and P>bo are also fitting parameters reported 
in Table II. Here again, a fourth order term has been 
included to achieve a better agreement in the fit of the 
energy which slightly deviates from the harmonic approx- 
imation at larger amplitude of polarizations. 



Abto -Bbto Cbto Abo Bbo 

-0.0239 9.369 10~ B -4.373 10" lu 0.122 1.514 10" 4 



TABLE II: Values of the parameters, used in Eq. (2) and Eq. 
(3), that result from a fit of the bulk internal energies under 
zero macroscopic field versus zero-field polarization shown in 
Fig. 1. Units required to produce energies in meV when the 
polarization enters in /iC/cm 2 . 



l^bto(-Pbto) — ^bto-Fbto + -Bbto-Pbto + Cbto-PbtO' 

(2) 

where AbtOj Bbto an d Cbto are fitting parameters, 
reported in Table II. The sixth order term was needed to 
ensure better agreement of the fit with the first-principles 
data. 




-60 -40 -20 20 40 60 
Polarization (|xC.cm~ 2 ) 



FIG. 1: (Color online) Energy as a function of spontaneous 
polarization for epitaxial bulks. The blue curve with squares 
corresponds to the tetragonal bulk BaO and the red curve 
with circles is the double well of tetragonal BaTiOs bulk. For 
each single point energy, the out-of-plane cell parameter was 
relaxed. 



B. Electrostatic energy cost 

In ferroelectric/incipient-ferroelectric super lattices, 
the incipient ferroelectric layer is expected to become 
polarized with nearly the same spontaneous polarization 
as the ferroelectric layer, resulting in a roughly homoge- 
neous polarization through the whole structure, 11 so that 
E e i ec in Eq. (1) will vanish. Here, in the BaTiOs/BaO 
super lattices, the BaO layer is less polarizable and is no 
more expected to have necessarily the same polarization 
than the BaTi03 layer. The polarization misfit gives 
rise to electric fields and produce an additional electro- 
static energy cost. Assuming that both layers arc homo- 
geneously polarized, its expression can be deduced from 
the expansion of the total energy of a dielectric with re- 
spect to the electric field £. 29 For the case of a single 
dielectric, 

E{£) = E{£ = 0) - £P° - ^e e°°£ 2 . (4) 

Generalizing Eq. (4) to the present context, where we 
have two dielectrics periodically repeated in the super- 
lattice, then the zero field energy can be written as 
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E(£ = 0) = mU BTO (P B > TO )+nU BO (P B > o ), (5) c{ ^ m) 

and the electrostatic energy for the superlattices can be 
written as 



-Eeiec(^BTO'^BO><?BTO,£Bo) — 

Q c pO ^BTO oo c2 
- "BTOC-BTO-fBTO ?, e e BTO^BTO 

^bo^bo-Pbo l^eoeBO^BO' ( 6 ) 



where SIbto an d f^BO are, respectively, the volumes of 
the BaTiC>3 and the BaO layers and can be estimated as 



Obto = tucbtoS, 
fi B o = nc BO S, 



(7a) 
(7b) 



with S the area of the surface of the unit cell [S = (3.84 
A) 2 ] and cbto an d cbo are approximated by c BTO and 

S) 30 

c bo- 

The electric fields inside the BaTi03 layer, £bto- an d 
the BaO layer, £bo, are n °t independent. From the 
short-circuit boundary conditions across the whole su- 
percell 



to cbto £ bto + n c B o £bo = 0, 



(8) 



the continuity of the normal component of the displace- 
ment field at the interface 



£ o£bto + -Pbto — £o£bo + -PbOi 
and the constitutive relations 

-Pbto = -Pbto + £ oXbto£bto> 
-Pbo = Pbo + £ oXbo^bo, 

we can arrive to expressions for £bto an d £bo ; 



£bto — 



n 



cbo (Pbto ~ -Pbo) 



(9) 



(10a) 
(10b) 



(11a) 



eo (n CBO £b°to + 171 c bto £bo) ' 
m cbto (-Pbto ~ p bo) 



eo (n cbo £b°to + m c bto ego) ' 

Replacing Eqs. (lla)-(llb) into Eq. (6), then the elec- 
trostatic energy reduces to 



E elec (n,m,P^ T0 ,P° ) = C(n,TO)(P° T0 -P° ) 2 , (12) 



where 



_S n m c B q cbto 

2e (to cbto ego + n c bo Cbto) 



(13) 



Here again, we will approximate cbto an d cbo by c BTO 
and Cg . 

Combining Eq. (5) and Eq. (12), the total energy for 
the electrostatic model for (BaTi03) m /(BaO)„ superlat- 
tices reads 



+ C(n,TO)(P° T0 -P° ) 2 . (14) 

C. Minimization of the total energy 

For each periodicity m/n, the ground state of the su- 
perlattice can be estimated by minimizing the total en- 
ergy Eq. (14) with respect to P BTO and P BO - Then, 
from the knowledge of the zero-field polarization, the 
constitutive relations [Eqs. (10a)-(10b)], and the elec- 
tric fields [Eqs. (lla)-(llb)], we can infer values for the 
total polarization inside each layer. In Table III we re- 
port the results of such minimization for a set of peri- 
odicities n/m. It is interesting to see that at a fixed 
number of BaO layers n, the model predicts a critical 
ratio m/n = 3.5 beyond which the system becomes fer- 
roelectric: 7/2, 14/4 and 21/6. For all of these ferro- 
electric states, the minimum of energy is reached when 
the polarization of the BaO layer is about 35% smaller 
than the polarization of the BaTi03 layer, giving rise to 
an electrostatic energy of the same order of magnitude 
than the internal energy. We notice that, considering 
the internal energy alone (i.e. assuming Pbto = ^bo an d 
E e i ec =0), we can also predict a critical ratio m/n ~ 5 
to get a ferroelectric state (11/2, 21/4 and 31/6). This 
critical ratio coming from the pure internal energies is 
however larger than the critical ratio (m/n ~3.5) cal- 
culated with the electrostatic energy, meaning that the 
quantity of BaTi03 to get ferroelectricity in the super- 
lattice is lowered by polarizing the BaO layer less than 
the BaTi03 layer, a situation distinct from that usually 
observed in ferroelectric/incipient-ferroelectric superlat- 
tices. As the polarizability of the insulating materials in- 
creases however, we will tend toward the situation where 

pO _ pO 

- r BTO~- r BO- 



IV. FIRST-PRINCIPLES RESULTS 

As a second step, we have performed DFT calcula- 
tions (i) to verify to which extent the predictions of the 
simple electrostatic model are verified, and (ii) to ana- 
lyze the eventual existence of instabilities in the phonon 
spectra, that would indicate the existence of structural 
phase transitions. 
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n 


m 


pU 

"bto 


pU 

■"BO 


U 


Eelec 


Etot 


2 


7 


3.2 


2.2 


-0.58 


0.57 


-0.01 


2 


8 


9.6 


6.3 


-6.69 


5.55 


-1.14 


2 


9 


13.0 


8.4 


-14.95 


11.03 


-3.92 


2 


10 


15.4 


9.9 


-24.77 


16.64 


-8.13 


2 


11 


17.3 


11.0 


-35.76 


22.17 


-13.59 


4 


14 


3.2 


2.2 


-1.17 


1.14 


-0.03 


4 


15 


7.2 


4.8 


-6.74 


6.03 


-0.71 


4 


16 


9.6 


6.3 


-13.38 


11.10 


-2.28 


4 


17 


11.4 


7.5 


-21.10 


16.43 


-4.67 


4 


18 


13.0 


8.4 


-29.90 


22.07 


-7.83 


6 


21 


3.2 


2.2 


-1.75 


1.71 


-0.04 


6 


22 


6.2 


4.1 


-7.09 


6.52 


-0.58 


6 


23 


8.1 


5.4 


-13.19 


11.48 


-1.71 


6 


24 


9.6 


6.3 


-20.07 


16.65 


-3.42 


6 


25 


10.9 


7.1 


-27.79 


22.11 


-5.68 



TABLE III: Polarizations P^to and -Pbo (mC/chi 2 ) mini- 
mizing the total energy of Eq. (14) for a range of different 
layer thicknesses n and m. The corresponding internal, elec- 
trostatic and total energies (respectively U, E e i ec and E to t) 
are also reported. Units of the energies in meV. 



A. Paraelectric reference 

In its paraelectric reference state with space group 
P4-/mmm (No 123), the in-plane atomic positions in each 
layer correspond to the bulk high symmetric positions. 
In the out-of-plane direction a mirror symmetry plane is 
present at the center of each layer. Since in the z direc- 
tion the atoms are not in their high symmetry positions, 
atomic relaxations take place and give rise to rumpling 
distortions. We recover the previous first-principles re- 
sults of Ref. 17, showing that the rumpling distortions 
are mainly localized at the interface between the layers 
and rapidly vanish at the center of thick enough layers, 
in good agreement also with the "locality principle" dis- 
cussed in Ref. 7. This interface effect in the superlattices 
involves also a modification of the interplanar distances 
between atomic planes with respect to those in the bulks, 
but this effect is also strongly located at the interface be- 
tween the two layers. 17 As shown in Table IV, this feature 
tends to compress the total out-of-plane cell parameter of 
the multilayer. In spite of the different periodicities con- 
sidered, a constant deviation is calculated between the 
bulk references and the relaxed superlattices confirming 
the fact that the differences with respect to the bulks are 
mainly located at the interface between the two layers. 

Since the ferroelectric instabilities are strongly sensi- 
tive to the volume and pressure, the reduction of the out- 
of-plane thicknesses in the superlattices could strongly 
influence the out-of-plane ferroelectric soft mode for the 
smallest periodicities. For example, in the 2/2 period- 
icity, this reduction with respect to the bulk reference 
equals 0.07 A, which corresponds to an equivalent pres- 
sure of about 1 GPa. 



n 


rn 




( "b ulk 


deviation 


2 


2 


13.43 


13.50 


-0.07 


2 


4 


21.45 


21.52 


-0.07 


2 


6 


29.46 


29.54 


-0.08 


2 


8 


37.48 


37.56 


-0.08 


4 


2 


18.91 


18.98 


-0.07 


4 


4 


26.93 


27.00 


-0.07 


4 


6 


34.95 


35.02 


-0.07 


4 


8 


42.95 


43.04 


-0.09 


4 


10 


50.99 


51.06 


-0.07 



TABLE IV: Out-of-plane cell parameter, c , for a range of dif- 
ferent layer thicknesses n and m. Cj olax represents the perpen- 
dicular cell parameter of the paraelectric relaxed case, c° u ik 
stands for the sum of the corresponding bulk unit cell, c° u ik 
= m Cg TO + § c bo- The last column reports the deviation 
between c° ulk and c^ olax . Units in A. 



B. Zone-center instability 

For two fixed BaO thicknesses, n=2 and n=4, the evo- 
lution of the square of the frequency of the lowest zone 
center TO mode along the out-of-plane direction (FE 2 ), 
w| , calculated in the relaxed paraelectric references, is 
shown in Fig. 2. This FE 2 mode, that after condensa- 
tion will be responsible of the ferroelectric transition, is 
strongly sensitive to the thickness of BaTi03 layer, as ex- 
pected since the ferroelectric instability is coming from 
BaTi03. As shown in Fig. 2, goes from high values 
at the smallest m and decreases linearly with the increase 
in the number of BaTiC>3 unit cells for both n=2 and 
n=4 BaO thicknesses. Beyond a critical m/n periodicity, 
7/2 and 9/4, the FE Z mode becomes unstable (negative 
w^ ). This means that, in spite of the non-ferroelectric 
character of BaO layer, the entire supercell can develop 
a ferroelectric instability when a critical ratio between 
m/n is reached. 




2 4 6 8 10 



m 

FIG. 2: (Color online) Square of the frequency of the low- 
est zone center out-of-plane transverse optic mode, uj^q, for 
n=2 (red, circles) and n=4 (blue, squares) as a function of 
the BaTiC>3 thickness, m. The phonon frequencies have been 
computed in the reference paraelectric structure. 
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This confirms the results predicted with the simple 
electrostatic model Eq. (14). The agreement is amaz- 
ingly good for n — 2 but the critical periodicities at which 
the ferroelectricity appears do not correspond to a con- 
stant ratio as inferred from the model: 7/2 and 9/4 from 
the first-principles phonon calculations instead of 7/2 and 
14/4 from the electrostatic model. Since the electrostatic 
model is built through the bulk soft mode eigendisplace- 
ments and the bulk Born effective charges (to estimate 
the polarization) , we now investigate to which extent this 
constitutes a reasonable approximation. 

To highlight the differences between the bulks and 
the superlattice we report in Table V the Born effective 
charges of individual atoms at the interface and in the 
middle of each layer of the superlattice in the paraelec- 
tric reference structure. In the middle of the BaO layer, 
the amplitude of the zz component of the Born effective 
charges are larger than their bulk values, a difference that 
is amplified when m increases and reduced when n in- 
creases. In the middle of BaTiC>3 the opposite behaviour 
is observed, with all the amplitudes of the Born effective 
charges smaller than in the BaTiC>3 bulk, and a trend to 
reach the bulk values when to increases while they are 
reduced when n increases. These global evolutions of the 
Born effective charges with n and m in the superlattices 
can be related to the natural disposition to recover the 
BaO or BaTiC>3 bulk values when respectively n or m 
increases. The evolution of the Born effective charges of 
Oi, Ti and O2/3 atoms at the interface are comparable to 
those in the middle of the BaTiC>3 layer, but with larger 
deviations from the bulk amplitudes. The Born effective 
charges of the Ba atom at the interface increases with m 
and decreases with n, but reach values for large to higher 
than the BaO or the BaTi03 bulks. The main conclusion 
that can be drawn from this discussion is that the Born 
effective charges of the atoms in the superlattice are sig- 
nificantly different than those at the bulk level, mainly 
for the smallest periodicities m/n. 

Since the Born effective charges of the atoms are mod- 
ified, we can expect a modification in the mode polarity 
of the superlattice FE Z instability with respect to the 
bulk. To compare the evolution of the ferroelectric mode 
for different periodicities m/n, we report in the Fig. 3(a) 
the dependence with respect the number of BaTi03 unit 
cells of the ratio between the oscillator strength and the 
volume of the unit cell (S/f2o) for n=2 and n=A. The 
oscillator strength of one mode is defined as 23 

S a p =^2Z*,aa'(VKa'T^2K',PP'VK'0', (15) 

KOc' k' P' 

where Z* are the individual atomic Born effective charge 
tensors and r\ is the phonon eigendisplacement vector. 

This S/Oo ratio increases strongly with m, evolving 
rapidly from values smaller than in bulk BaO at m=2 
to larger values. At the considered thicknesses, S/f2o is 
much smaller than for bulk BaTi03 but we can expect 
that it will reach the BaTi03 bulk value at large enough 
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-1.97 
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TABLE V: zz components of the Born effective charges (e) of 
atoms in the reference paraelectric BaTiOs /BaO superlattices 
with different periodicities m/n. The first two lines refer the 
Ba and O atoms in the middle of the BaO layer. Lines 3rd to 
6th refer to the Ba, Ti and O atoms at the interface. Lines 
7th to 10th refer to the Ba, Ti and O atoms in the middle 
of the BaTiOs layer. Oi refers to the O in the BaO plane 
while O2/3 stands for the two equivalent oxygens in the Ti02 
planes. 



to. It is clear here that the modification of the S/f2o 
ratio comes either from a modification of the Born effec- 
tive charge and/or from a modification of the eigenvector 
n. As shown in Table V, the Born effective charges only 
tends to increase smoothly with to and cannot be re- 
sponsible to the strong modification of S/flg shown in 
Fig. 3(a). 

In Fig. 3(b) we report the evolution of the square of 
the norm of the phonon eigendisplacement vector rf 1 with 
respect to to. Here, rj 2 increases strongly with to, which 
can be related to a strong enhancement of the ferroelec- 
tric eigenvector r\ with the BaTiOs thickness. Therefore, 
according to Eq. (15), the evolution of the S/^o ratio 
can be mainly attributed to the modification of the fer- 
roelectric eigenvector when to is modified. 



In summary, this analysis points out that the sim- 
ple electrostatic model based on the bulk quantities can 
fail to reproduce the first-principles results because of 
a strong modification of the ferroelectric eigenvector in 
the superlattice for the smallest periodicities. This fer- 
roelectric mode is strongly sensitive to the thickness of 
the BaTiOa layer in the superlattice and, surprisingly 
does not reach the BaTiOs bulk ferroelectric soft mode, 
even for the largest m/n considered. Compared to other 
previously studied superlattices like SrTi03/BaTi03,the 
modification of the soft-mode pattern could be ampli- 
fied here due to the fact that the BaTiOs and the BaO 
layer are not sharing the same perovskite structure. This 
teaches us that the quantitative agreement of the model 
for to = 2 could be at least partly fortuitous and that 
model prediction have to be considered with caution. 
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FIG. 3: (Color online) Evolution of the S/f2o ratio (a) and 
T] 2 = J]] VkVk (b) of the ferroelectric soft mode for n=2 (red, 
circles) and n=4 (blue, squares) as a function of the BaTiC>3 
thickness m. The horizontal dashed lines correspond to the 
BaTi03 and BaO bulk limit S/fio values which are respec- 
tively equal to 43 1CF 7 and 5 10 -7 a.u. 



C. Phonon dispersion curves 




The phonon dispersion curves between V (0,0,0), X 
(|,0,0) and M (|,|,0) points were also calculated, and 
the results are shown in Fig. 4 for periodicities 2/2, 4/2 
and 6/2. At the smallest thickness (2/2), no instability 
is present, neither at T, X nor M point. This confirms 
that the ground state of the 2/2 multilayer is the non- 
polar PJ^/mmm phase. However, for the 4/2 periodicity 
an instability appears at the X zone-boundary point with 
a frequency of 42i cm , and its magnitude is amplified 
for larger to, with a calculated frequency of 92i cm" 1 for 
6/2. Moreover, in Fig. 4(b) a low frequency mode is also 
observed at the M zone boundary point for the 4/2 thick- 
ness. This mode becomes unstable for the 6/2 periodic- 
ity, but with an amplitude smaller than the instability at 
the X point. The branch responsible for the M point in- 
stability is the same than the one responsible for X point 
instability and its dispersion is flat between the X and 
M points while it shifts rapidly to a positive frequency 
when going from X or M point to the T point. Similar 
behaviors are also observed for the to/4 periodicities. An 
X point instability is also observed for m >4 and its am- 
plitude is also amplified with m, with esentially the same 
imaginary frequencies as in the n = 2 case (42i cm -1 
for 4/4 and 91i cm -1 for 6/4). In conclusion, for both 
thicknesses of the BaO layer, zone boundary instabilities 
appear before the zone center ferroelectric becomes soft. 



FIG. 4: (Color on line) Phonon dispersion curves for 
m/n=2/2 (a), 4/2 (b), and 6/2 (c). The color of the different 
branches are assigned according to the contribution of each 
chemical species to the dynamical matrix eigenvector (red for 
the Ba atom, green for the Ti atom, and blue for the O atom) . 



To clarify the atomic motions related to the X point 
instability, we display a schematic picture of the eigendis- 
placement pattern in Fig. 5 for the 6/4 periodicity. Since 
the instability is located at the X point, atoms in consec- 
utive cells along the [100] direction move out-of-phase. 
As we can see on Fig. 5, these eigendisplacements can be 
decomposed into two components: (i) polar distortions 
along the z direction, only for the atoms along the Ti- 
O chains, and (ii) in-plane motions, only for Ba and O 
atoms which are along the Ba-0 chains parallel to the 
Ti-0 chains. The amplitudes of these in-plane displace- 
ments are however much smaller than the out-of-plane 
displacements, meaning that the Ti-0 distortions along 
the z direction dominate the total motions (blue-green 
color of the unstable branch at the X point in Fig. 4). 
Moreover, on the two BaO atomic planes corresponding 
to the two mirror planes of symmetry in the z direction, 
the in-plane motions completely disappear. Additionally, 
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we can notice that the polar distortions along the z di- 
rection remain in the BaO layer but only for the atoms 
aligned with the Ti-0 chain, in order to preserve the flow 
of large polarization currents along those chains. 

Since in consecutive unit cells along the [100] direc- 
tion the polar displacements are in opposite directions, 
these X point instabilities are related to an antiferro- 
electric distortion. To highlight these issues, we report 
also in Fig. 5(b) schematic vectors showing the direction 
of main local dipolar moments. From this picture, it is 
clear that along the z directions we have antiferroelectric 
distortions, with alternating chains of up and down po- 
larization. For the in-plane Ba and O motions, the atoms 
move only in the [100] direction, giving rise to dipolar mo- 
ments oriented along the x direction. Between the two 
mirror symmetry planes, these in-plane dipolar moments 
keep the same orientation, but it is inverted after cross- 
ing these symmetry planes. Interestingly, this inversion 
does not take place at the interface of the supercell but 
in the middle of each layer. 
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FIG. 5: (Color online) (a) Atomic displacements of the X 
point unstable mode for the 6/4 thickness. Two consecutive 
cells along the [100] direction are shown. Atoms are repre- 
sented by balls (O in red, Ti in green and Ba in blue). Dashed 
lines represent the position of the mirror symmetry planes of 
the paraelectric reference structure. Dotted lines are the po- 
sition of the interfaces, (b) Schematic representation of the 
direction of polarization resulting from local polar distortions. 



V. DISCUSSION 



The thickness evolution of the phonon dispersion 
curves of the BaTiOs/BaO superlattice, as discussed in 
Sec. IV C, can be understood in relationship with those 
of bulk BaTiC>3. At the bulk level, BaTiC>3 was shown to 
exhibit a chain-like ferroelectric instability in real space, 
with weak inter-chain interactions, yielding ferroelectric 
modes similarly unstable at V, X and M points. 31 ' 32 The 
ferroelectric mode at T is highly polar. In the super- 
lattice, due to the difficulty to polarize BaO, it will be 
associated to a huge depolarizing field which will strongly 
harden it. At the opposite, the zone-boundary ferroelec- 
tric modes at X and M, are globally non-polar. They will 
therefore not induce any macroscopic depolarizing field 
and will so keep a much stronger tendency to be unsta- 
ble than the T mode. Consequently, the zone-boundary 
modes will always be more unstable than the T mode 
in the superlattice, yielding an anti-ferroelectric ground- 
state. It is only in the limit of thick BaTiC>3 layers, for 
which the depolarizing field tends to vanish, that we will 
recover the usual ferroelectric ground-state. It is also 
worth to notice that, for n = 2 and 4, the instability at 
X similarly appears at m = 4 , which corresponds to what 
was previously reported as the approximative length of 
the correlation volume (to ss 4 — 5) required to induce a 
ferroelectric instability. 31 ' 33 

The evolution of the ground state of the superlattice 
with the ticknesses of the layers can therefore be summa- 
rized as follows: (i) For the smallest BaTi03 thicknesses 
to < 4, the system does not exhibit any ferroelectric in- 
stability and the supercell remains paraelectric. (ii) For 
larger BaTiC>3 thicknesses (to > 4) but small ratio m/n, 
the ferroelectric instability at T is suppressed by the ef- 
fect of the depolarizing field; there is only an instability 
at X and the system is antiferroelectric. (iii) As the ra- 
tio m/n increases, the ferroelectricity will progressively 
expand from X to T and we can expect the formation of 
ferroelectric domains of increasing sizes, (iv) Only when 
to >> n, the superlattice will tend to the purely ferro- 
electric ground-state. 

We notice that the present BaTiOs/BaO super- 
lattices present similarities with Ruddlesden-Popper 
A n+ iB n 0^ n+ i materials which can be viewed as the 
stacking of alternating perovskite and rocksalt layers. 18 
In the latter, however, there is no continuity of the Ti-0 
chains from one perovskite block to the next one. This 
will prevent the possibility of huge polarization currents 
along the stacking direction as those associated to the 
feroelectric mode in BaTiOs/BaO system and could ex- 
plain why Ruddlesden-Popper materials do not develop 
any tendency to become ferroelectric with polarization 
aligned along the stacking direction. 19 
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VI. CONCLUSIONS 

In this paper we have studied, from first-principles, su- 
perlattices combining m unit cells of BaTi03 with n unit 
cells of BaO. Although BaTi03 is ferroelectric at the bulk 
level, we have shown that the superlattice can present an 
antiferroelectric ground-state and explained that by the 
hardening of the ferroelectric mode with respect to the 
antiferroelectric mode due to depolarizing field issues. 

This behavior is quite general. While epitaxial strain 
was already previously reported to be able to modify the 
competition between ferroelectric and antiferrodistortive 
instabilities, we propose that playing with the depolariz- 
ing field in ferroelectric/insulator superlattices is another 
practical way to tune the competition between polar and 
non-polar instabilites. We have demonstrated here the 
possibility to engineer an antiferroelectric ground-state 
in a system combining ferroelectric and antiferroelectric 



instabilities. We anticipate that in systems with com- 
peting ferroelectric and antiferrodistortive instabilities, it 
could similarly favor an antiferrodistortive ground-state. 
We hope that these results will motivate further studies 
of feroelectric/insulator superlattices.. 
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